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Higgs boson production with subsequent decay to photons was searched for in 

the data collected by the DELPHI detector at centre-of-mass energies between 
183 GeV and 209 GeV, corresponding to a total integrated luminosity of nearly 
650 pb _1 . No evidence for a signal was found, and limits were set on h°Z° and 
h°A° production with h° decay to photons. These results were used to exclude 
regions in the parameter space of fermiophobic scenarios of Two Higgs Doublet 
Models. 



(Accepted by Euro. Phys. J. C) 



J.Abdallah 25 , P.Abreu 22 , W.Adam 51 , P.Adzic 11 , T.Albrecht 17 , T.Alderweireld 2 , R.Alemany-Fernandez 8 
T.Allmendinger 17 , P.P.Allport 23 , U.Amaldi 29 , N.Amapane 45 , S.Amato 48 , E.Anashkin 36 , A.Andreazza 28 , S.Andringa 22 
N.Anjos 22 , P.Antilogus 25 , W-D.Apel 17 , Y.Arnoud 14 , S.Ask 26 , B.Asman 44 , J.E.Augustin 25 , A.Augustinus 8 , P.Baillon 8 
A.Ballestrero 46 , P.Bambadc 20 , R.Barbier 27 , D.Bardin 16 , G.J.Barker 17 , A.Baroncelli 39 , M.Battaglia 8 , M.Baubillier 25 
K-H. Becks 53 , M.Begalli 6 , A.Behrmann 53 , E.Ben-Haim 20 , N.Benekos 32 , A.Benvenuti 5 , C.Berat 14 , M.Berggren 25 
L.Berntzon 44 , D.Bertrand 2 , M.Besancon 40 , N.Bcsson 40 , D.Bloch 9 , M.Blom 31 , M.Bluj 52 , M.Bonesini 29 , M.Boonekamp 40 
P.S.L. Booth 23 , G.Borisov 21 , O.Botner 49 , B. Bouquet 20 , T.J.V.Bowcock 23 , I.Boyko 16 , M.Bracko 43 , R.Brenner 49 
E.Brodet 35 , P.Bruckman 18 , J.M.Brunet 7 , L.Bugge 33 , P.Buschmann 53 , M.Calvi 29 , T.Camporesi 8 , V.Canale 38 , F.Carena 8 
N.Castro 22 , F.Cavallo 5 , M.Chapkin 42 , Ph.Charpentier 8 , P.Checchia 36 , R.Chierici 8 , P.Chliapnikov 42 , J.Chudoba 8 
S.U.Chung 8 , K.Cieslik 18 , P.Collins 8 , R.Contri 13 , G.Cosmc 20 , F.Cossutti 47 , M.J.Costa 50 , D.Crcnnell 37 , J.Cuevas 34 
J.D'Hondt 2 , J.Dalmau 44 , T.da Silva 48 , W.Da Silva 25 , G. Delia Ricca 47 , A.Dc Angelis 47 , W.De Boer 17 , C.De Clercq 2 
B.Dc Lotto 47 , N.Dc Maria 45 , A.Dc Min 36 , L.de Paula 48 , L.Di Ciaccio 38 , A.Di Simonc 39 , K.Doroba 52 , J. Drees 53 ' 8 
M.Dris 32 , G.Eigen 4 , T.Ekelof 49 , M.Ellert 49 , M.Elsing 8 , M.C.Espirito Santo 22 , G.Fanourakis 11 , D.Fassouliotis 11 ' 3 
M.Feindt 17 , J.Fernandez 41 , A.Ferrer 50 , F.Ferro 13 , U.Flagmeyer 53 , H.Foeth 8 , E.Fokitis 32 , F.Fulda-Quenzer 20 , J.Fuster 50 
M.Gandclman 48 , C.Garcia 50 , Ph.Gavillet 8 , E.Gazis 32 , R.Gokieli 8 ' 52 , B.Golob 43 , G.Gomcz-Ccballos 41 , P.Goncalvcs 22 
E.Graziani 39 , G.Grosdidicr 20 , K.Grzclak 52 , J.Guy 37 , C.Haag 17 , A.Hallgren 49 , K.Hamacher 53 , K.Hamilton 35 , S.Haug 33 
F. Hauler 17 , V.Hedberg 26 , M.Hennecke 17 , H.Herr 8 , J.Hoflfman 52 , S-O. Holmgren 44 , P.J.Holt 8 , M.A.Houlden 23 
K.Hultqvist 44 , J.N.Jackson 23 , G.Jarlskog 26 , P.Jarry 40 , D.Jeans 35 , E.K.Johansson 44 , P.D.Johansson 44 , P.Jonsson 27 
C.Joram 8 , L.Jungermann 17 , F.Kapusta 25 , S.Katsancvas 27 , E.Katsoufis 32 , G. Kernel 43 , B.P.Kersevan 8 ' 43 , U.Kerzel 17 
A.Kiiskinen 15 , B.T.King 23 , N.J.Kjaer 8 , P.Kluit 31 , P.Kokkinias 11 , C.Kourkoumelis 3 , O.Kouznetsov 16 , Z.Krumstein 16 
M.Kucharczyk 18 , J.Lamsa 1 , G.Lcdcr 51 , F.Lcdroit 14 , L.Leinonen 44 , R.Leitner 30 , J.Lemonne 2 , V.Lepeltier 20 , T.Lesiak 18 
W.Licbig 53 , D.Liko 51 , A.Lipniacka 44 , J. H. Lopes 48 , J.M.Lopez 34 , D.Loukas 11 , P.Lutz 40 , L.Lyons 35 , J.MacNaughton 51 
A.Malek 53 , S.Maltezos 32 , F.Mandl 51 , J.Marco 41 , R.Marco 41 , B.Marechal 48 , M.Margoni 36 , J-C. Marin 8 , C.Mariotti 8 
A.Markou 11 , C.Martinez-Rivero 41 , J.Masik 12 , N.Mastroyiannopoulos 11 , F.Matorras 41 , C.Matteuzzi 29 , F.Mazzucato 36 
M.Mazzucato 36 , R.Mc Nulty 23 , C.Meroni 28 , E.Migliore 45 , W.Mitaroflf 51 , U.Mjocrnmark 26 , T.Moa 44 , M.Moch 17 
K.Moenig 8 ' 10 , R.Monge 13 , J.Montenegro 31 , D.Moraes 48 , S.Moreno 22 , P.Morettini 13 , U.Mueller 53 , K.Muenich 53 
M.Mulders 31 , L.Mundim 6 , W.Murray 37 , B.Muryn 19 , G.Myatt 35 , T.Myklebust 33 , M.Nassiakou 11 , F.Navarria 5 
K.Nawrocki 52 , R.Nicolaidou 40 , M.Nikolenko 16 ' 9 , A.Oblakowska-Mucha 19 , V.Obraztsov 42 , A.Olshcvski 16 , A.Onofrc 22 
R.Orava 15 , K.Osterberg 15 , A.Ouraou 40 , A.Oyanguren 50 , M.Paganoni 29 , S.Paiano 5 , J.P.Palacios 23 , H.Palka 18 
Th.D.Papadopoulou 32 , L.Pape 8 , C.Parkes 24 , F.Parodi 13 , U.Parzefall 8 , A. Passed 39 , O.Passon 53 , L.Peralta 22 
V.Perepelitsa 50 , A.Pcrrotta 5 , A.Petrolini 13 , J.Piedra 41 , L.Pieri 39 , F.Pierre 40 , M.Pimenta 22 , E.Piotto 8 , T.Podobnik 43 
V.Poireau 8 , M.E.Pol 6 , G.Polok 18 , V.Pozdniakov 16 , NPukhacva 2 ' 16 , A.Pullia 29 , J.Ramcs 12 , A.Read 33 , P.Rebccchi 8 
J.Rehn 17 , D.Reid 31 , R.Rcinhardt 53 , P.Rcnton 35 , F.Richard 20 , J.Ridky 12 , M.Rivero 41 , D.Rodriguez 41 , A.Romero 45 
P.Ronchcsc 36 , P.Roudeau 20 , T.Rovelli 5 , V.Ruhlmann-Kleidcr 40 , D.Ryabtchikov 42 , A.Sadovsky 16 , L. Salmi 15 , J. Salt 50 
C. Sander 17 , A.Savoy-Navarro 25 , U.Schwickerath 8 , A.Segar 35 , R.Sekulin 37 , M.Siebel 53 , A.Sisakian 16 , G.Smadja 27 
O.Smirnova 26 , A.Sokolov 42 , A.Sopczak 21 , R.Sosnowski 52 , T.Spassov 8 , M.Stanitzki 17 , A.Stocchi 20 , J.Strauss 51 , B.Stugu 4 
M.Szczekowski 52 , M.Szeptycka 52 , T.Szumlak 19 , T.Tabarelli 29 , A.C.TafTard 23 , F.Tegenfeldt 49 , J.Timmermans 31 
L.Tkatchcv 16 , M.Tobin 23 , S.Todorovova 12 , B.Tomc 22 , A.Tonazzo 29 , P.Tortosa 50 , P.Travnicek 12 , D.Treille 8 , G.Tristram 7 
M.Trochimczuk 52 , C.Troncon 28 , M-L.Turluer 40 , I.A.Tyapkin 16 , P.Tyapkin 16 , S.Tzamarias 11 , V.Uvarov 42 , G.Valenti 5 
P. Van Dam 31 , J. Van Eldik 8 , A. Van Lysebetten 2 , N.van Remortel 2 , I. Van Vulpcn 8 , G.Vegni 28 , F.Veloso 22 , W.Venus 37 
P.Verdier 27 , V.Verzi 38 , D.Vilanova 40 , L.Vitalc 47 , V.Vrba 12 , H.Wahlen 53 , A.J.Washbrook 23 , C.Weiser 17 , D.Wicke 8 
J.Wickcns 2 , G.Wilkinson 35 , M.Winter 9 , M.Witek 18 , O.Yushchenko 42 , A.Zalewska 18 , P.Zalewski 52 , D.Zavrtanik 43 
V.Zhuravlov 16 , N.I.Zimin 16 , A.Zintchenko 16 , M.Zupan 11 



1 Department of Physics and Astronomy, Iowa State University, Ames IA 50011-3160, USA 
2 Physics Department, Universiteit Antwerpen, Universiteitsplein 1, B-2610 Antwerpen, Belgium 
and IIHE, ULB-VUB, Pleinlaan 2, B-1050 Brussels, Belgium 

and Faculte des Sciences, Univ. de l'Etat Mons, Av. Maistriau 19, B-7000 Mons, Belgium 
3 Physics Laboratory, University of Athens, Solonos Str. 104, GR-10680 Athens, Greece 
4 Department of Physics, University of Bergen, Allegaten 55, NO-5007 Bergen, Norway 
5 Dipartimcnto di Fisica, Universita di Bologna and INFN, Via Irncrio 46, IT-40126 Bologna, Italy 
6 Centro Brasileiro de Pcsquisas Fisicas, rua Xavier Sigaud 150, BR-22290 Rio de Janeiro, Brazil 
and Depto. de Fisica, Pont. Univ. Catolica, C.P. 38071 BR-22453 Rio de Janeiro, Brazil 
and Inst, de Fisica, Univ. Estadual do Rio de Janeiro, rua Sao Francisco Xavier 524, Rio de Janeiro, Brazil 
7 College de France, Lab. de Physique Corpusculaire, IN2P3-CNRS, FR-75231 Paris Cedex 05, France 
8 CERN, CH-1211 Geneva 23, Switzerland 

9 Institut de Recherches Subatomiques, IN2P3 - CNRS/ULP - BP20, FR-67037 Strasbourg Ccdcx, France 
10 Now at DESY-Zeuthen, Platanenallee 6, D-15735 Zeuthen, Germany 

"Institute of Nuclear Physics, N.C.S.R. Dcmokritos, P.O. Box 60228, GR-15310 Athens, Greece 

12 FZU, Inst, of Phys. of the C.A.S. High Energy Physics Division, Na Slovancc 2, CZ-180 40, Praha 8, Czech Republic 
13 Dipartimento di Fisica, Universita di Genova and INFN, Via Dodccancso 33, IT-16146 Gcnova, Italy 
14 Institut des Sciences Nucleaires, IN2P3-CNRS, Universite de Grenoble 1, FR-38026 Grenoble Cedex, France 
15 Helsinki Institute of Physics, P.O. Box 64, FIN-00014 University of Helsinki, Finland 

16 Joint Institute for Nuclear Research, Dubna, Head Post Office, P.O. Box 79, RU-101 000 Moscow, Russian Federation 
17 Institut fur Expcrimcntellc Kcrnphysik, Universitat Karlsruhe, Postfach 6980, DE-76128 Karlsruhe, Germany 
18 Institute of Nuclear Physics PAN,U1. Radzikowskiego 152, PL-31142 Krakow, Poland 

19 Faculty of Physics and Nuclear Techniques, University of Mining and Metallurgy, PL-30055 Krakow, Poland 
20 Universite de Paris-Sud, Lab. de l'Accelerateur Lineaire, IN2P3-CNRS, Bat. 200, FR-91405 Orsay Cedex, France 
21 School of Physics and Chemistry, University of Lancaster, Lancaster LAI 4YB, UK 
22 LIP, 1ST, FCUL - Av. Elias Garcia, 14-1°, PT-1000 Lisboa Codex, Portugal 
23 Department of Physics, University of Liverpool, P.O. Box 147, Liverpool L69 3BX, UK 
24 Dept. of Physics and Astronomy, Kelvin Building, University of Glasgow, Glasgow G12 8QQ 

25 LPNHE, IN2P3-CNRS, Univ. Paris VI ct VII, Tour 33 (RdC), 4 place Jussieu, FR-75252 Paris Ccdcx 05, France 
26 Department of Physics, University of Lund, Solvegatan 14, SE-223 63 Lund, Sweden 
27 Univcrsitc Claude Bernard dc Lyon, IPNL, IN2P3-CNRS, FR-69622 Villeurbanne Cedex, France 
28 Dipartimento di Fisica, Universita di Milano and INFN-MILANO, Via Ccloria 16, IT-20133 Milan, Italy 
29 Dipartimcnto di Fisica, Univ. di Milano-Bicocca and INFN-MILANO, Piazza della Scicnza 2, IT-20126 Milan, Italy 
30 IPNP of MFFj Charles Univ., Areal MFF, V Holesovickach 2, CZ-180 00, Praha 8, Czech Republic 
31 NIKHEF, Postbus 41882, NL-1009 DB Amsterdam, The Netherlands 

32 National Technical University, Physics Department, Zografou Campus, GR-15773 Athens, Greece 
33 Physics Department, University of Oslo, Blindern, NO-0316 Oslo, Norway 
34 Dpto. Fisica, Univ. Oviedo, Avda. Calvo Sotelo s/n, ES-33007 Ovicdo, Spain 
35 Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK 
36 Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, IT-35131 Padua, Italy 
37 Rutherford Appleton Laboratory, Chilton, Didcot OX11 OQX, UK 

38 Dipartimento di Fisica, Universita di Roma II and INFN, Tor Vergata, IT-00173 Rome, Italy 
39 Dipartimcnto di Fisica, Universita di Roma III and INFN, Via della Vasca Navale 84, IT-00146 Rome, Italy 
40 DAPNIA/Service de Physique des Particules, CEA-Saclay, FR-91191 Gif-sur-Yvette Cedex, France 
41 Institute de Fisica de Cantabria (CSIC-UC), Avda. los Castros s/n, ES-39006 Santander, Spain 
42 Inst. for High Energy Physics, Serpukov P.O. Box 35, Protvino, (Moscow Region), Russian Federation 
43 J. Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia and Laboratory for Astroparticle Physics, 
Nova Gorica Polytechnic, Kostanjcviska 16a, SI-5000 Nova Gorica, Slovenia, 
and Department of Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
44 Fysikum, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

45 Dipartimento di Fisica Sperimentale, Universita di Torino and INFN, Via P. Giuria 1, IT-10125 Turin, Italy 
INFN,Sezione di Torino, and Dipartimento di Fisica Teorica, Universita di Torino, Via P. Giuria 1, 
IT-10125 Turin, Italy 
47 Dipartimento di Fisica, Universita di Trieste and INFN, Via A. Valerio 2, IT-34127 Trieste, Italy 

and Istituto di Fisica, Universita di Udine, IT-33100 Udine, Italy 
48 Univ. Federal do Rio de Janeiro, C.P. 68528 Cidade Univ., Ilha do Fundao BR-21945-970 Rio dc Janeiro, Brazil 
49 Department of Radiation Sciences, University of Uppsala, P.O. Box 535, SE-751 21 Uppsala, Sweden 
50 IFIC, Valcncia-CSIC, and D.F.A.M.N., U. de Valencia, Avda. Dr. Moliner 50, ES-46100 Burjassot (Valencia), Spain 
51 Institut fur Hochcncrgicphysik, Osterr. Akad. d. Wissensch., Nikolsdorfcrgasse 18, AT-1050 Vienna, Austria 
52 Inst. Nuclear Studies and University of Warsaw, Ul. Hoza 69, PL-00681 Warsaw, Poland 
53 Fachbcrcich Physik, University of Wuppertal, Postfach 100 127, DE-42097 Wuppertal, Germany 



1 Introduction 

In the Standard Model (SM), the decay of the Higgs boson to photons is mediated 
by heavy charged particles (namely W ± bosons and top-quarks). The corresponding 
branching ratio BR(h° — > 77) is below 0.1%. At LEP, this search is thus motivated by 
extensions of the SM. Many of the proposed models may enhance the BR(h° — > 77), 
either by enlarging the /i°77 coupling [1] or by reducing the coupling of the Higgs boson 
to fermions [2]. 

We investigate the fermiophobic scenario of Two Higgs Doublets Models (2HDM) [2], in 
which the lightest CP-even Higgs boson decays to photons and can be produced together 
with a Z° or a CP-odd Higgs boson, A . The Higgs decay modes analysed in this paper are 
h° — > 77 and A — > bb or A — > h°Z°. However, we consider also results for h° — »■ A°A° 
and for long-lived A , as described in [3]. All LEP2 data with centre-of-mass energies 
above 180 GeV are analysed and the results reported here update those in [3]. 

The results on h°Z° production with h° — > 77 have been interpreted in several frame- 
works: previous analyses of DELPHI data can be found in [4] and results from the other 
LEP experiments can be found in [5]. In all cases, h°Z° production with subsequent 
h° — > 77 decay has been used as a benchmark in the search for particles with SM Higgs- 
like couplings to bosons but no coupling to fermions. 

In sections 2 and 3, the 2HDM fermiophobic scenario and the data samples used will 
be introduced; the general selection of events with isolated photons will then be presented 
(in section 4), before going into the analysis dedicated to each of the h° boson production 
mechanisms and the different final states. The results obtained for h°Z° production and 
h°A° production (in sections 5 and 6, respectively) are then combined to exclude regions 
in the parameter space of the 2HDM fermiophobic scenario, in section 7. 

2 The 2HDM fermiophobic scenario 

General Two Higgs Doublets Models without explicit CP violation are characterized by 
five physical Higgs bosons: two neutral CP-even bosons (h° and H°) , two charged bosons 
{H^) and one neutral CP-odd boson (A ). Together with the masses, the important 
parameters describing 2HDMs are the mixing angle in the neutral CP-even sector (a) 
and the ratio of the vacuum expectation values of the two Higgs doublets (tan fS) . 

The couplings of the Higgs doublets to fermions could be realized in different ways, 
one possibility is that only one of the doublets couples to fermions. The coupling of the 
lightest CP-even boson to a fermion pair is then proportional to cos a. If a = tt/2, this 
coupling vanishes and h° becomes a fermiophobic Higgs [2]: it decays to pairs of other 
Higgs bosons or massive gauge bosons when kinematically allowed, or to two photons in 
a large region of the parameter space. 

In 2HDMs, the main mechanisms for the production of neutral Higgs bosons at LEP 
are e + e~ — ► h°Z° and e + e~ — ► h°A°, both proceeding via Z° exchange. The two processes 
have complementary cross-sections proportional to sin 2 5 and cos 2 5, respectively, where 
5 = a — [3. The sin S factor rescales the h°Z°Z vertex with respect to the SM one. 

It must be noticed that there are two different Higgs potentials which conserve CP [2] 
(referred to as potential A and B), each of them a function of seven parameters - the 
masses of the Higgs bosons, the a and (3 angles and the sum of the squares of the vacuum 
expectation values for the two doublets. The choice of the potential does not affect 
the Higgs interactions with gauge bosons or fermions but leads to different Higgs-Higgs 



couplings and thus different phenomenologies. In particular, the decay width of h° — > 77, 
to which diagrams with H^ loops significantly contribute, can be changed. 

For potential A, the branching ratio of the lightest CP-even Higgs to two photons 
BR(/i° — > 77) depends on M h o, and weakly on sin 2 <5. Under the assumption that 
Mho ~ 1 TeV/c 2 and that h° does not decay to other Higgs particles, the BR(/i° — » 77) 
for sin 2 5=1 can be obtained using the SM couplings to bosons and no coupling to 
fermions. For other values of sin 2 5, it increases with respect to this benchmark model. 
For potential B, the BR(/i° — > 77) depends also on M^o and M H ± and there can be large 
cancellations between the several loop contributions for some values of these parameters. 
In both cases, h° — > A A is dominant when M^o > 2M^o. 

The widths for the A tree level decays (to //, Z°h° and W^H^) are independent 
of the potential but depend on M^o, M^o and sin 2 S. For low values of sin 2 5, also A 
becomes fermiophobic and light A bosons can thus be long lived. 

The low values of sin 2 5 lead also to different mass configurations, according to the 
choice of potential. For potential A they imply M h o ~ and for potential B M h o ~ M^o. 

3 Data samples 

The data analysed were taken by the DELPHI detector at LEP in the years from 
1997 to 2000. The corresponding average centre-of-mass energies (y/s) and integrated 
luminosities (C) are shown in Table 1. 



year 


1997 


1998 


1999 


2000 


VS (GeV) 
C (pb- 1 ) 


182.6 
49.3 


188.6 
153.0 


191.6 
25.1 


195.5 
76.0 


199.6 

82.7 


201.6 
40.2 


205.0 
80.0 


206.5 
59.2 


206.8 
81.8 



Table 1: Average centre-of-mass energies and integrated luminosities of the analysed 
samples. The sample with ^=206.5 GeV corresponds to the data collected after the 
damage on the TPC. 



In the year 2000, the centre-of-mass energies ranged from 200 GeV to 209 GeV, while 
most of the luminosity was collected at around 205 GeV and 207 GeV. In the last part of 
the running, DELPHI suffered irreparable damage to one of the sectors of the main track- 
ing device, the Time Projection Chamber (TPC), representing 1/12 of the acceptance. 
These data were analysed separately, to isolate any systematic difference. 

The DELPHI apparatus and performance are described in detail in [6,7]. The tracking 
system of DELPHI consisted of the TPC and a Vertex Detector (VD) closest to the beam 
pipe and was supplemented by extra tracking detectors, the Inner and Outer Detectors 
in the barrel region, and two Forward Chambers. It was embedded in a magnetic field of 
1.2 T, aligned parallel to the beam axis. The most important subdetectors for the present 
analysis were the electromagnetic calorimeters covering different polar angle regions 1 : the 
luminosity monitor (STIC) for 6 below 11°, the Forward ElectroMagnetic Calorimeter 
(FEMC) between 11° and 35° and the High density Projection Chamber (HPC) above 
42°. The regions in between the FEMC and the HPC and the HPC intermodular division 
at 6 — 90° were equipped with hermeticity counters - scintillators covered with lead, so 
that photons could also be tagged there. In the azimuthal intermodular divisions of the 

The polar angle, 8, is defined in relation to the beam axis. In all cases the complementary value (180° — 9) is also 
assumed. The azimuthal angle, <j>, is the angle in the plane perpendicular to the beam direction. 



HPC, at mod(4>, 15°)=7.5°, the detection of photons could be complemented by the use of 
the Hadronic CALorimeter (HCAL). The hadronic calorimeter covered 98% of the total 
solid angle, down to 11°, and the whole detector was surrounded by muon drift chambers. 
The major hardware change with respect to the description in [7] was the inclusion of 
the Very Forward Tracker [8] which extended the coverage of the Vertex Detector down 
to a polar angle of 11°. Together with new tracking algorithms, and new alignment 
and calibration procedures, this led to an improved track reconstruction efficiency in the 
forward regions of DELPHI. The tracking algorithms for the barrel part of DELPHI were 
also changed to recuperate efficiency in the damaged TPC sector. 

Events corresponding to the SM processes and Higgs production signals were fully 
simulated for each data set. The main background processes were generated with KK2f [9] 
(for e + e~ — ► qqi'j)), KoralZ [10] (for e + e~ — ► vv(pj) and e + e~ — ► l + l^('~f)) and Bhwide [11] 
(for the Bhabha scattering), all of which include a detailed description of the initial 
state radiation. These processes constitute irreducible sources of background, and thus 
their accurate description is a crucial point in the analysis. In a previous paper [3], the 
QQ(l) background had been simulated with Pythia 6.1 [12]; KK2f has a more accurate 
description of the initial state photon radiation and gives a 30% increase in the cross- 
section for the production of qq with two photons at high polar angle. It also has a better 
match to the DELPHI data as shown in the next sections. 

Other SM processes give smaller contributions to the studied channels. All the 
4-fermion processes (including neutral and charged currents) were generated as described 
in [13], using the WPHACT generator [14], and complemented with samples of multipe- 
ripheral production of e + e~ f f events generated with BDK/BDKRC [15] in the regions 
dominated by virtual photon collisions. The QED e + e~ — > 77 process and Compton 
events were generated with RADCOR [16] and TEEG [17], respectively. The hadroniza- 
tion and fragmentation in hadronic final states was performed by Pythia [12]. 

Signal processes corresponding to the different analysed final states were generated 
according to Pythia 6.1 [12] at the different average centre-of-mass energies, for different 
values of the Higgs bosons masses (from 3 to 120 GeV/c 2 for the h° boson and 10 to 
170 GeV/c 2 fro the A boson). They were cross-checked with a dedicated generator 
HZHA [18]. 

All the generated data sets at the different centre-of-mass energies were passed through 
the DELPHI simulation and the same reconstruction chain as data [7]. Dedicated samples 
of all the above processes were used to simulate the effect of the damaged TPC sector. 
Its impact on the analyses was found to be negligible within the statistical uncertainties. 

4 Selection of events with isolated photons 

The selection of events with isolated photons is common to all analyses and is described 
in this section. First a general selection was applied and then isolated photons and leptons 
were reconstructed. 

Only events with visible energy in the polar angle region above 20° greater than 0.1 ^/s 
were accepted. In addition, all events were required to contain at least one charged or 
neutral object with energy above 5 GeV in the same polar angle region. This vetoed 
most of the contamination from virtual photon collision events. 

Charged particles were classified as 'good' if they had measured momentum greater 
than 0.1 GeV/c and their tracks were reconstructed with impact parameters to the in- 
teraction point below 4 cm in the transverse plane and 4/sin 9 cm in the beam direction. 



Energy deposits in the calorimeters unassociated to charged particles were required to be 
above 0.3 GeV to be classified as neutral particles. 

The reconstruction of photons was done in several steps, starting from the showers in 
the electromagnetic calorimeters. The procedure described in [7] was followed to iden- 
tify tight photon candidates in the HPC. This algorithm selects electromagnetic energy 
deposits with shower profiles compatible with those of photons. However, showers close 
to the HPC divisions were accepted as loose photon candidates even if they failed the 
transverse shower profile criteria (and/or the longitudinal shower profile criteria if their 
energy was above 25 GeV). In the forward region, all STIC energy deposits with polar 
angle satisfying 9 < 11° were taken to be tight photon candidates 2 . An algorithm was 
used to correct for the effects of photon conversion and shower development in the detec- 
tor material in front of the FEMC, as explained below. Electromagnetic deposits close 
in space were clustered together and the association with reconstructed charged parti- 
cle tracks was used for electron/photon discrimination. Tight photon candidates were 
required to have no association to VD track elements, nor to signals from different com- 
binations of other tracking detectors (depending on the shower polar angle). Care was 
taken to exclude those tracks which were likely to come from the development of showers 
outside the calorimeter. Loose photon candidates were allowed to have two associated 
tracks. In addition, the ratio of electromagnetic energy to the total energy around the 
cluster - in an angular region defined by |A0| < 15° and |A0| < min(15°, 6° cot 9 c i ust er) - 
was required to be above 90%. 

The identification of isolated photons started from the candidates defined above and 
used a double cone centred around the photon axis, as explained below. Only isolated 
photons with energies above 5 GeV were considered. The selection criteria depended on 
the topology of the event and were as follows. 

For the topology with photons only, the total energy inside a cone with half-angle of 
10° was associated to the photon, while, to ensure isolation, the total energy between 10° 
and 12° was not permitted to exceed 5 GeV. 

For the other topologies, the criteria were different for loose and tight photon candi- 
dates. For loose photon candidates, the half-angle of the inner cone was of 5° and the 
energy between 5° and 15° was not permitted to exceed 1 GeV. In the case of photons 
tightly identified by the shower profile analysis, no further association was done and only 
the external cone was kept. The angle of the external cone, a, was varied according to 
the energy of the photon candidate, down to 3° for E 1 > 90 GeV. The energy limit inside 
the cone was rescaled to sin a/ sin 15° GeV, but one energetic particle was accepted inside 
this cone (and excluded from the energy calculation). 

The identification of isolated photons in the barrel region required also that there was 
no HPC layer with more than 90% of their electromagnetic energy, while the hadronic 
energy depositions above 3 GeV must be concentrated in the first layer of the HCAL. 

The reconstruction of isolated leptons followed the same double cone criteria, starting 
from good charged particles with momentum above 4 GeV/c. The total charged particle 
momentum and energy deposition inside the inner cone around the charged particle track 
were associated, and an external cone used to ensure isolation. Isolated charged particles 
associated to signals in the muon chambers, and for which the ratio between the energy 
deposited in calorimeters and the measured momentum was less than 20%, were identified 
as tight muon candidates. Electromagnetic showers associated to good charged particles, 
reconstructed with algorithms similar to those of photon reconstruction, were considered 
as tight or loose electron candidates. 

2 Energy depositions below 3° were discarded from the events, to avoid contamination from off-momentum beam 
electrons. 



The separation between electrons and photons converting in the tracking system relied 
on the information from the VD and its association to isolated charged particles. A VD 
track element was defined as at least two signals in different layers of the detector, which 
were associated to an isolated particle if aligned within 3° of its azimuthal direction (10° in 
the forward region). Isolated charged particle tracks not associated to VD track elements 
were considered as candidates for photons converting in the tracking system. At most 
one converted photon candidate was allowed per event, except for the / + /~77 topologies 
where no recovery of converted photons was performed. 

5 Search for h°Z° production 

5.1 Two photons and jets 

In h°Z° production, the final state with highest branching fraction is the one corre- 
sponding to the Z° decay to hadrons, and h° decay to photons. The preselection of this 
topology required that at least six good charged particles were present and the visible 
energy in the polar angle region above 20° was greater than 0.2-y/s. All selected charged 
and neutral particles, except for isolated photons and leptons, were clustered into two 
jets using the DURHAM jet algorithm [19]. Two isolated photons were required in the 
event, and extra isolated particles were allowed only if their transverse momentum with 
respect to one of the jets was less than 20 GeV/c. 

The main irreducible background to this search was qq production with two photons 
coming predominantly from initial state radiation: at LEP2, about half of the qq events 
were radiative return events, with an effective collision energy around M z o. Very small 
contributions from other sources, namely W + W~ production, were also present. 

After the selection of the two isolated photons in the event, further requirements on 
their isolation and polar angle reduced both Initial and Final State Radiation (ISR and 
FSR) contributions. The most energetic photon had to have at least 15% of the beam 
energy, the minimum transverse momentum of the photons with respect to the jets had 
to be above 7.5 GeV/c and the minimum polar angle above 30°. 

To improve the resolution on the invariant mass of the photon pair, a kinematic fit 
imposing energy-momentum conservation and constraining the jet-jet invariant mass to 
the Z° mass was performed according to [20] . The main inputs to the fit are the directions 
and energies of the photons, which are considerably less well measured if the photons are 
not well contained in the calorimeters, oe, the relative error on the measured energy 
EZf as , was obtained from calorimeter resolution studies. It was set to 100% for pho- 
tons reconstructed in the calorimeter boundaries, and events in which both photons had 
oe = 100% were rejected. For the other events, a partial x 2 was constructed as: 

2 , / Tpjit _ rnmeas 

A'w - 1 - \ *> / 'yi ' 



Y" — f ~ 1± — (I) 



ndof 2 .f^ ^ a E E lt 

If both photons were well contained in the calorimeters, it was required that xL,/ndo/ 
was below 5. If one of the photons was not well contained, the information from the jets 
had a larger weight and it was the global x 2 /ndof of the fit (defined in [20]) that was 
required to be below 5. 

At this stage, the contribution from non-gg(7) events becomes almost negligible. The 
characteristics of the radiative return events were then used to reduce the background 
further. In most of these events, the photons are in the forward part of the detector and, 



usually, one of the photons carries most of the energy necessary to bring the Z° boson 
on-shell: 

E ret = ^— ^. (2) 

Only events in which the energy difference between the two photons was lower than 
0. 70 E ret were kept. Figure f shows the sum and difference of the two photon energies 
normalized to the radiative return energy. In the final sample, one of the photons was 
required to be inside the HPC acceptance. 

Table 2 shows, for the data collected in the year 2000, the evolution with the selec- 
tion criteria of the number of selected data events, corresponding SM background and 
efficiencies for two Higgs signals of different masses. The contribution of the irreducible 
qq background is shown separately. The numbers of selected events in all data and simu- 
lated samples at the final selection level are shown in Table 5. The reconstructed jet-jet 
invariant masses and the fitted 77 masses are also shown in Figure 1. 

5.2 Two photons and two charged leptons 

The production of h°Z° with subsequent decays h° — ► 77 and Z° — ► l + l~ was analysed 
for all the three charged lepton flavours. A common preselection required the presence of 
two isolated photons and two isolated leptons in events with at most five good charged 
particles and visible energy in the polar angle region above 20° greater than 0.2-^/s. The 
most energetic photon was required to have energy greater than 0.1^/s and the polar angle 
of the leptons was required to be above 20° (except for the ones identified as muons). 

The background processes for these channels include e + e~ — > Z /'-/* — > l + l~ and the 
t-channel Bhabha scattering, which has a very large cross-section. To reduce the number 
of events with FSR, the minimum transverse momentum of each photon with respect 
to any of the leptons was required to be greater than 5 GeV/c. The contribution from 
events with ISR, namely from radiative return to the Z°, was reduced requiring that both 
photons had polar angles above 30°, that there was at least one photon in the HPC, and 
imposing that the difference between the energies of the two photons was below 0.7E ret 
(see equation 2). 

A kinematic fit imposing energy-momentum conservation and using the measured 
directions of the four particles in the event was performed, to rescale the energies of 
leptons and photons. Partial x 2 /ndof associated to the leptons or to the photons were 
defined, as in equation 1, for the photon energies and lepton momenta; events were 
accepted only if one of the two was below 10. Only events for which the fitted di-lepton 
invariant mass was between 60 and 120 GeV/c 2 were kept. 

Table 3 presents the number of selected events in the year 2000 data and the corre- 
sponding SM background expectations, at each selection level, together with the efficien- 
cies for two Higgs masses. The selection efficiencies for the different leptonic decays of 
the Z° vary between 13% for t + t~ and 38% for fi + fi~ for a Higgs mass of 100 GeV/c 2 . 
The final invariant mass spectra are shown in Figure 2, and the corresponding numbers 
of selected events in each data set are given in Table 5. 

5.3 Two photons and missing energy 

h Z production with the Z° decay into neutrino pairs leads to purely photonic final 
states with large missing energy. All purely photonic candidate events were allowed to 
have at most 5 good charged particles, none of them associated to VD track elements or to 



signals in the muon chambers. Charged particles not associated to energy deposits above 
5 GeV had to have momenta below 5 GeV/c and the minimum transverse momentum 
of each individual particle with respect to the other particles reconstructed in the same 
hemisphere had to be greater than 5 GeV/c. 

Cosmic rays leaving energy deposits in the calorimeters are an important source of 
background for final states with photons and missing energy. Most cosmic ray events 
crossing the tracking system were removed by requiring all charged particles in the event 
to satisfy the impact parameter selection criteria defined in section 4. Cosmic ray events 
crossing the detector outside the tracking devices were vetoed by requiring that the total 
unassociated energy deposition in the event was less than 10% of the energy of the isolated 
photons. Events in which more than 98% of the energy of the photon candidates was 
reconstructed from depositions in the hadronic calorimeters were also vetoed. In addition, 
the direction of photons reconstructed in the HPC was required to be consistent within 
25° with the hypothesis that the particle was coming from the primary vertex. 

After this preselection, only events with two photons reconstructed with polar angles 
above 25° and well contained either in the FEMC or in the HPC were kept. 

The main background process for this channel is the production of vv pairs (either 
through Z° or W ± boson exchange) with the emission of ISR photons. Even though the 
cross-section is low, these events constitute an irreducible source of background. The 
background from the QED reaction e + e~ — > 77(7), which has a large cross-section, 
was much reduced by appropriate kinematic requirements. Most of the QED events were 
vetoed by imposing that the angle between the two photons was lower than 178° and that 
the energy of the most energetic photon was below 0.4-^/i. Additional criteria were applied 
to events with photons in the less efficient regions of the electromagnetic calorimeters. If 
one of the photons was either within 1.5° of an azimuthal modular division of the HPC, 
or its polar angle corresponded to the HPC edges, (defined by the intervals [42°, 44°] and 
[87°, 88°]), the acoplanarity 3 was required to be above 3° and 5°, respectively. If there 
was one converted photon in the FEMC, the acoplanarity was also required to be above 
3°, in order to account for the larger deflection of charged particles in the forward region. 
Any signal observed in the hermeticity counters in between the FEMC and the HPC was 
required to be within 20° of a reconstructed photon. 

A kinematic fit was then performed imposing the Z° mass on the invisible system, 
after requiring that the missing mass of the selected event was larger than 20 GeV/c 2 . 
The global y 2 /ndof [20] resulting from the fit was required to be below 5. 

The results of the selection in the year 2000 data set are summarized in Table 4, where 
the SM expectations were corrected for trigger efficiencies 4 . The invariant mass spectra at 
the final selection level are shown in Figure 2, and the corresponding numbers of selected 
events in each data set are given in Table 5. 

5.4 Results on h°Z° production 

A good agreement between data and SM expectations was found in all analysed chan- 
nels, with a total of 54 events selected and 52±1 expected. The numbers of selected 
events per channel and centre-of-mass energy are shown in Table 5. 

Limits on the production cross-section of h°Z° with h° — > 77, as a function of M h o, 
were obtained by combining all the channels and centre-of-mass energies using the Mod- 

3 The acoplanarity between two objects is defined as the complement of the angle between them in the plane transverse 
to the beam direction. 

The trigger efficiencies for events with two photons in the barrel part of DELPHI were of 96%-98%; in the forward 
part, as for other final states, they were of >~100%. In the other channels the effect of this correction is thus negligible. 



ified Frequentist Likelihood Ratio method [21], taking into account the measured and 
expected M 77 invariant mass distributions. The charged lepton flavours were added in a 
single channel for each centre-of-mass energy. 

The signal samples generated with Pythia 6.1 were cross-checked with samples gen- 
erated with HZHA and found to be compatible within the statistical uncertainty of 1%. 
Different fragmentation/hadronization models were also compared, and found to have 
negligible impact in the qq'j'j selection. Extra systematic effects could come from the 
parameterization of the expected invariant 77 mass distributions. They were found to be 
small when compared to the statistical uncertainty on the total background expectations. 
To account for all effects, a systematic error of ± 3% was assigned to the signal efficiency 
used in the limit calculation. 

Figure 3 shows, as a function of the Higgs boson mass, the 95% Confidence Level (CL) 
upper limit on BR(h° — > 77) times the ratio of the h°Z° production cross-section to the 
SM one. This ratio is equal to sin 2 5 in 2HDMs (see section 2). For the range of masses 
studied, and taking the model in [2] as a reference, the width of the Higgs boson is always 
below the mass resolution of the analyses 5 . 

In a model where the Higgs couplings to bosons have SM values, but the couplings to 
fermions vanish, a 95% CL lower limit on the h° boson mass is given by the intersection 
of the cross-section limit and the prediction for BR(h° — > 77) (also shown in Figure 3), 
at 104.1 GeV/c 2 . The expected limit in this case is 104.6 GeV/c 2 . 

6 Search for h°A° production 

6.1 Two photons and two b-jets 

The search for h°A° — > 7766 is very similar to the one for h°Z° — > 7755 (see section 
5.1). Similar criteria were applied, differing mainly in the kinematic fit performed. In 
the h°A° case less constraints were used: the jet-jet invariant mass was left free and an 
additional ISR photon was allowed in the beam direction. In fact, if M^o + M^o < M^o, 
a significant part of the signal events may correspond to radiative returns to the Z° 
with subsequent Z° — ► h°A°. Since M40 has to be fitted without constraints, the global 
X 2 /ndof becomes an important selection criterion: it was required to be below 10 in all 
cases. Tighter requirements on partial xi^/ndof and global x 2 / n d°f were still apllied as 
discussed for the h°Z° case. 

In the next stage, it was required that the energies of the two photons satisfied 
l-Eyi — E j2 \ < 0.70|i? 7 i + E l2 \- Again, at least one of the photons had to be in the HPC 
and the combined b-tagging of the event, as defined in [22], was required to be above 
-2. Figure 4 shows the agreement of the b-tagging variable in data and SM simulation, 
together with the mass distributions obtained for the bb and 77 pairs. 

The total numbers of selected data and background events in the analysis of the year 
2000 data are shown in Table 6. The non-gg background is larger compared to the case 
with the Z° mass constraint: it consists mainly of semi-leptonic W + W~ events. Table 6 
shows also the efficiency obtained for e + e~ — » h°A° with M A o = M z o and M h o = 100 
GeV/c 2 . The numbers of selected events at the final analysis level are given in Table 8. 

5 For M h a = 100 GeV/c 2 , the mass resolutions were of 1% for the leptonic channel, 2.5% in the hadronic channel, and 
3.5% in the missing energy channel. 



6.2 Three or more photons and jets 

Within the framework of general 2HDM, the main decay channel of A is bb. However, 
when kinematically allowed, the decay to h°Z° becomes dominant in a large region of 
parameter space. This can lead to final states with a fermion pair and four photons. Only 
the Z° decay to quark pairs, which has the highest branching ratio, was analysed. At 
least three isolated photons were required and at most one converted photon was allowed. 

Two of the photons had to be above 30° in polar angle and have a transverse momen- 
tum with respect to any of the jets greater than 7.5 GeV/c. The jet-jet invariant mass 
was required to be in the range 50 to 130 GeV/c 2 . The difference between the energy 
of any of the photons and the energy carried by all the others had to be below 0.33E ret 
(defined in equation 2). Finally, one of the photons was required to be in the HPC, and 
at least three photons were required to be above 15° in polar angle. 

The numbers of selected events in the 2000 data sample and the corresponding SM 
expectations are given in Table 7. The efficiencies for the signal were as high as 60% for 
high masses, and above 40% for most of the considered range. Since only three photons 
were required, no mass reconstruction was attempted and the result is based on the 
numbers of selected events at this stage (given in Table 8). 

6.3 Results on h°A° production 

The numbers of selected events in each data set and for both channels are shown in 
Table 8. Limits were extracted with the same algorithm as for h°Z°. In the btryy channel 
the reconstructed mass values were used, while for the 99777 (7) there was no attempt 
to reconstruct the masses, and a pure counting experiment was performed. Searches in 
both channels had similar sensitivity to the signal. 

To extract upper limits on cos 2 8 (defined in section 2), reference cross-sections com- 
puted with HZHA [18] were used. For the h° — > 77 decay, the branching ratio computed 
with sin 2 5 = 1 and shown in Figure 3 was used. This represents a conservative assump- 
tion since BR(/i° — > 77) increases for lower sin 2 S. Depending on the point in parameter 
space under consideration, the limit was derived assuming a 100% branching fraction for 
the dominant decay channel of the A (bb or h°Z°). The limits are almost independent 
of the other model parameters, if the masses M#o and Mjj± are high enough not to open 
new decay channels as A — »■ H^W^. 

Figure 5 shows the limits on sin 2 5 as a function of M^o, obtained for two different 
Mjsp values - they strongly depend on the mass hypothesis used for A . 

7 Exclusion in the 2HDM parameter space 

The combination of the results on h°Z° and h°A° production is illustrated in Figure 5. 
The upper limits on sin 2 5 for a given M^o and the upper limits on cos 2 5 for a given 
(M^o ,Mjyo) pair are combined to exclude regions in the (M^M^o) plane for all S values. 

Figure 6 shows the 95% CL exclusion in the plane (M h o,M^o) divided into separate 
regions according to the different kinematics: 

I. M40 > M h o + M z o 

The main decay channels in this region are h° — > 77 and A — ► h°Z°. Thus the 77Z 
and 99777(7) topologies are the relevant ones. The vertex A°h°Z° is proportional to cos 5 



10 



and thus the decay A — > h°Z° is dominant only for sin 2 5 < 0.95. However, for higher 
sin 2 S values, these mass combinations are excluded by the limits on h°Z° production. 

II. M h o/2 < M A o < M h o + M z o 

h° — > 77 and A — > bb are the dominant decays in most of the parameter space. 
However, the decay A — »■ bb is suppressed as sin 2 5 approaches 0, and the A becomes 
stable, and thus invisible, for sin 2 5 < 10" 6 . The exclusion in this region comes from the 
search for h°A° production, the results of a 77 + E mis analysis (described in [3]) being 
used. This allows all the mass region excluded by the combination of the h°Z° and 6677 
analyses for higher sin 2 S values also to be excluded for low sin 2 5. 

III. 2M b < M A o < M h o/2 

The decay h° — > A A , kinematically allowed, becomes dominant. The A still decays 
to bb, as long as 2Mb < M A o, and this gives rise to 6-fermion final states. This decay is 
relevant also for non-fermiophobic 2HDMs, and the results published by DELPHI in [23] 
are used. For very low sin 2 5 there is a region of totally invisible final states for which 
the results on the invisible Z° width, obtained at LEP1, are used, as explained in [3]. 

The LEP1 results on the total Z° width are also used to cover the region where the 
two masses are below the reach of the present analysis [3] . 

These exclusions are valid for all the allowed parameter space in the case of potential 
A, but only for M#+ >500 GeV/c 2 or sin 2 5 > 0.1 in the case of potential B. For lower 
values of both Mh+ and sin 2 S, the BR(/i° — > 77) can vanish due to cancellation of the 
different loop contributions. 

For the two potentials considered, the exclusions obtained for invisible A (which are 
more restrictive than the ones for visible final states, in region III) apply only to small 
mass regions. This is because sin 2 S < 10 -6 implies M h o ~ in potential A and a band of 
M h o ~ M40 in potential B. Both mass bands are in principle outside region III but, since 
their precise width is not known, the region of exclusion for invisible A is also shown. 

8 Conclusions 

Around 650 pb _1 of LEP2 data collected by DELPHI, at centre-of-mass energies be- 
tween 183 and 209 GeV, were analysed in the search for Higgs bosons decaying into 
photons. In the context of 2HDM, both h°Z° and h°A° production were searched for, 
and a large variety of final states involving photons and fermions was considered. No 
evidence for new physics was found. 

Lower limits were set on the mass of a particle with Higgs-like couplings to bosons 
and decaying to two photons. In a model where the Higgs couplings to bosons have 
SM values, but the couplings to fermions vanish, a 95% CL lower limit on the h° boson 
mass is set at 104.1 GeV/c 2 ; the expected limit in this case is 104.6 GeV/c 2 . Exclusions 
at 95% CL were also derived in the mass plane (M h o,M A o) of the fermiophobic 2HDM 
scenario. 
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Selection 


Data 


Background 


qq 


eff 50 


effioo 


2 photons 


338 


341.5 ±9.5 


296.5 ±9.2 


60% 


60% 


no-ISR/FSR 


23 


22.9 ±0.3 


21.2 ±0.3 


39% 


45% 


Zfit 


12 


12.4 ±0.2 


12.3 ±0.2 


31% 


41% 


AE 


9 


8.7 ±0.2 


8.6 ±0.2 


26% 


40% 


HPC 


9 


8.4 ±0.2 


8.3 ±0.2 


26% 


40% 



Table 2: h°Z° — > 'yyqq selection in the year 2000: the evolution of the data sample with 
the analysis cuts is shown for the full 2000 data set, compared to the total expected 
background (and the qq background) and the efficiencies for Higgs masses of 50 GeV/c 2 
and 100 GeV/c 2 . The errors on the backgrounds are statistical only. The absolute 
statistical error on the efficiencies is 1-2%. 



Selection 


Data 


Background 


eff 50 


effioo 


2 photons 


24 


28.8 ±1.1 


36% 


36% 


no-ISR/FSR 


5 


5.3 ±0.5 


22% 


31% 


E-p fit ± M z o 





1.6 ±0.2 


17% 


27% 



Table 3: h°Z° — > 77/+/^ selection in the year 2000: the evolution of the data sample 
with the analysis cuts is shown together with the total expected background (with the 
corresponding statistical errors) and the efficiencies for Higgs masses of 50 GeV/c 2 and 
100 GeV/c 2 (with 1-2% statistical error). 



Selection 


Data 


Background 


QED 


z/z/77 


eff 5 o 


effioo 


preselection 


1077 


1121.2 ±7.9 


1110.6 ±7.9 


10.6 ±0.8 


60% 


64% 


2 7 s : 9 y i [35°, 42°] 


879 


924.5 ±7.2 


915.7 ±7.1 


8.8 ±0.7 


50% 


56% 


a 77 < 178° 


206 


200.0 ±3.3 


191.2 ±3.2 


8.8 ±0.7 


50% 


56% 


QED veto 


13 


9.0 ±0.7 


1.2 ±0.3 


7.8 ±0.7 


48% 


53% 


Zfit 


6 


6.6 ±0.6 


0.5 ±0.2 


6.1 ±0.6 


44% 


50% 



Table 4: h°Z° — > 77Z/Z/ selection in the year 2000: the evolution of the data sample 
with the analysis cuts is shown, together with the total expected background (and the 
QED and z/z/77 expectations) and corresponding statistical errors. The efficiencies for 
Higgs masses of 50 GeV/c 2 and 100 GeV/c 2 are shown in the last two columns, and have 
associated statistical errors of 1-2%. 
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3 
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0.29 ±0.07 


4 


1.58 ±0.22 


207 


3 


3.11 ±0.14 





0.69 ±0.17 





2.58 ±0.31 


TOT 


32 


27.15 ±0.47 


3 


5.03 ±0.55 


19 


19.49 ±1.06 



Table 5: Number of selected events for h°Z° production and corresponding background 
expectations for the three topologies and all data samples considered. The errors on the 
background are statistical only. 



Selection 


Data 


Background 


qq 


effioo 


no-ISR/FSR 


23 


22.9±0.3 


21.2±0.3 


53% 


fit 


19 


20.1±0.3 


19.1±0.3 


49% 


AE 


14 


17.1±0.3 


16.1±0.2 


49% 


HPC 


14 


16.5±0.3 


15.5±0.3 


49% 


b-tag 


7 


6.9±0.2 


6.4±0.2 


46% 



Table 6: h°A° — > 7766 analysis: the comparison of data and MC events selected at each 
analysis level for the year 2000 data. The 55(77) contribution is shown separately. The 
efficiency for the selection of a h°A° signal with M^o = 100 GeV/c 2 and M40 = M^o 
produced at y / i=206 GeV is shown in the last column (the corresponding statistical 
error being of 1-2%). 



Selection 


Data 


Background 


qq 


eff 


no-ISR/FSR 


6 


4.22 ±0.14 


4.00 ±0.14 


66% 


Z° mass 


3 


3.31 ±0.13 


3.07 ±0.12 


66% 


AE 


2 


1.85 ±0.10 


1.65 ±0.08 


64% 


HPC 


1 


0.86 ±0.07 


0.73 ±0.06 


62% 



Table 7: h°A° — ► 55777(7) analysis: the comparison of selected data and MC at each 
analysis level for the year 2000 data is shown. The efficiency for the selection of a h°A° 
signal with M h o = 35 GeV/c 2 and M^o = 135 GeV/c 2 produced at ^=206 GeV is shown 
in the last column (the corresponding statistical error being of 1-2%). 
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2.43 ±0.09 


1 


0.32 ±0.04 


206.5 


1 


1.85 ±0.09 





0.20 ±0.03 
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2.59 ±0.10 





0.34 ±0.05 


TOT 


21 


22.76 ±0.42 


3 


2.75 ±0.15 



Table 8: Selected events for h°A° production and corresponding background expectations 
for the two topologies and all data samples considered. The errors on the background 
are statistical only. 
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Figure 1: /i°Z° — > jjqq. The top plots show the sum (left) and difference (right) of the 
two photon energies divided by the radiative return energy before the cut on the energy 
difference (indicated by the arrow). The bottom plots show the reconstructed invariant 
jet-jet mass (left) and the fitted Higgs mass (right). The full data set (dots) is compared 
to the SM background (shaded area) and two h°Z° signals: with 71^0=50 GeV/c 2 and 
M fe o = 100 GeV/c 2 (thick lines), shown with arbitrary normalization. 
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Figure 2: h°Z° — > 77/+/^ (top) and /i°Z° — > ^vv (bottom). The invariant masses M 77 
at the last selection level for all analysed data (dots) are compared to the SM background 
expectations (shaded histogram), and h°Z° signals of 50 and 100 GeV/c 2 (top), shown 
with arbitrary normalization. 
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Figure 3: 95% CL upper limit on the h°Z° production cross-section x BR(/i° — > 77) 
normalized to the SM value. Both the observed (full line) and the expected limits (dashed 
line) are shown. Also shown is the fermiophobic BR(h° — ► 77) (dotted line), obtained 
by keeping the SM couplings of the Higgs to boson pairs and setting the h°ff couplings 
to 0. A 95% CL mass limit is given by the intersection of the cross-section limit and the 
BR(h° -> 77) curve at 104.1 GeV/c 2 (104.6 GeV/c 2 expected). 
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Figure 4: /i°y4° — > 7766 analysis. The reconstructed A and /i° masses at the final 
selection level are shown in the top plots for data and expected SM background in all 
data sets. The lower plots show the fitted jet-jet mass for events selected in both the 
h°A° and h°Z° analyses and the distribution of the b-tag variable just before the cut at 
-2. 
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Figure 5: 95% CL limits on the h°Z° and h°A° production cross-sections expressed in 
terms of sin 2 5. Values above the h°Z° curve and below the h°A° curves (shown for two 
different M40 values: 50 GeV/c 2 and 115 GeV/c 2 ) are excluded. 
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Figure 6: The shaded areas correspond to regions excluded at 95% CL for all 5 values. 
The plot is divided into regions according to the dominant decay modes of h° and A , as 
explained in the text. The exclusions from LEP 1 data, based on the total and invisible 
width of the Z°, are also shown. 



